Phase diagram of a model for a binary mixture of nematic molecules on a Bethe lattice 
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We investigate the phase diagram of a discrete version of the Maier-Saupe model with the inclusion of addi- 
tional degrees of freedom to mimic a distribution of rodlike and disklike molecules. Solutions of this problem 
on a Bethe lattice come from the analysis of the fixed points of a set of nonlinear recursion relations. Besides the 
fixed points associated with isotropic and uniaxial nematic structures, there is also a fixed point associated with a 
biaxial nematic structure. Due to the existence of large overlaps of the stability regions, we resorted to a scheme 
to calculate the free energy of these structures deep in the interior of a large Cayley tree. Both thermodynamic 
and dynamic-stability analyses rule out the presence of a biaxial phase, in qualitative agreement with previous 
mean-field results. 



I. INTRODUCTION 

The recent characterization of biaxial nematic phases in 
thermotropic liquid crystals [1] has renewed the theoretical in- 
terest in mechanisms leading to macroscopic biaxiality in such 
systems. In general, two possibilities exist [|2|]: (i) Either the 
systems are composed of intrinsically biaxial mesogens or (ii) 
the systems consist of a mixture of rodlike and disklike meso- 
gens. Although unquestionable experimental evidence of bi- 
axiality has only surfaced for systems of the first type, there 
have been hints from both experiments Ot] and computer sim- 
ulations [4] that mixtures can escape the curse of phase seg- 
regation and form stable biaxial nematic phases under the ap- 
propriate conditions. In this paper we introduce a framework 
which will allow future theoretical investigations of what such 
conditions might be, starting from microscopic models and 
obtaining analytical results which go beyond mean-field cal- 
culations. As a first application, we show that a very sim- 
ple choice of the interaction potential between mesogens in a 
binary mixture leads to a biaxial state which is unstable to- 
wards segregation, in agreement with previous mean-field re- 
sults @]. 

The inclusion of additional degrees of freedom in the 
Maier-Saupe model has been used to mimic a mixture of rod- 
like and disklike molecules, and to provide an explanation 
for the appearance of a biaxial nematic structure olla]. Ac- 
cording to some investigations for a mean-field Maier-Saupe 
model restricted to a discrete set of orientations, which we 
call the Maier-Saupe-Zwanzig (MSZ) model, the inclusion 
of a fixed distribution of shape variables leads to a stable bi- 
axial nematic phase. The biaxial region of the phase diagram 
is separated by critical lines from two distinct uniaxial ne- 
matic phases [6], in qualitative contact with some experimen- 
tal phase diagrams for a lyotropic liquid mixture This 
quenched polymorphism, however, which is generally used in 
solid-state systems, may not be adequate for liquids and liq- 
uid crystalline systems, with relatively short relaxation times, 
which may be better represented by thermalized degrees of 
freedom. We then carried out a mean-field investigation of 
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the analogous MSZ model, with an annealed (or thermalized) 
distribution of shape variables [8]. At the mean-field level, we 
have shown that a biaxial solution of this thermalized problem 
is still present, but it becomes thermodynamically unstable, 
and therefore physically unacceptable. 

In this article we report an analysis of a similar MSZ model 
on the deep interior of a Cayley tree, also known as a Bethe 
lattice. The main purpose of this investigation is the analysis 
of the global phase diagrams under the effects of fluctuations 
and short-range correlations neglected in the simple mean- 
field picture, but allowed for by the structure of the Cayley 
tree. In Section 2, we define the Zwanzig or discrete version 
of the Maier-Saupe model. The statistical problem is formu- 
lated in terms of a set of nonlinear recursion relations, whose 
fixed points correspond to solutions deep in the interior of a 
large tree 1191 1 1 Oil - This simple MSZ model displays a first- 
order transition between a disordered and a uniaxial ordered 
nematic phase, which is indicated by an overlap of the tem- 
perature ranges of (dynamic) stability of the fixed points as- 
sociated with disordered and ordered structures. The location 
of this first-order boundary comes from the application of an 
ingenious scheme, due to Gujrati Hill , which leads to the cor- 
rect thermodynamic free energy corresponding to each attrac- 
tor, and avoids the well-known pathologies associated with the 
surface of the Cayley tree IU2I1 . In the infinite-coordination 
limit, we recover the well-known mean-field results. In Sec- 
tion 3, we introduce the MSZ model for a thermalized binary 
mixture of rodlike and disklike molecules. Again, the prob- 
lem is formulated as a set of recursion relations, whose fixed 
points correspond to the physical solutions. Besides the at- 
tractors associated with the disordered and two ordered uni- 
axial structures, we find an additional fixed point associated 
with a biaxial nematic structure. There is a large region of 
overlap of (dynamic) stability of the two nematic uniaxial at- 
tractors. However, we find that the attractor associated with 
biaxial nematic phase is always dynamicaly unstable. Also, 
by using Gujrati's method, we show that the fixed point asso- 
ciated with the biaxial structure corresponds to a larger free 
energy with respect to the coexisting uniaxial attractors. This 
thermodynamic analysis and the (dynamic) stability analysis 
rule out the physical presence of an equilibrium biaxial phase 
in the Bethe lattice, in qualitative agreement with the previous 
mean-field results. 



II. THE MSZ MODEL ON A CAYLEY TREE 



The Maier-Saupe model is described by the Hamiltonian 



;tf = -A£ £ sfsf, 

(ij) f i,v=x,y,z 



(1) 



where A is a positive parameter, the first sum is over neigh- 
boring sites on a lattice of N sites, and Sf is an element of the 
traceless tensor 



sf = i(3<«r-s, 



(2) 



where is the /j-component of the unit vector associated with 
the orientation of a molecular aggregate on site i. In the Maier- 
Saupe-Zwanzig (MSZ) model we restrict the unit directors to 
the Cartesian axes, «,■ = ±(1,0,0), ±(0,1,0), ±(0,0,1). The 
canonical partition function is given by 



Z = £ exp 



{»,•} 



(3) 



where p is the inverse of the temperature. Since S^ v is invari- 
ant under the transformation n, = —Hi, the problem is reduced 
to the analysis of a three-state model, which can be shown to 
lead to the same qualitative features of the continuous Maier- 
Saupe model 11311 . 

We now formulate the MSZ model on a Cayley tree. In Fig- 
ure [T] we represent some generations (or layers) of a (rooted) 
Cayley tree of coordination q = 3. Also, we illustrate one of 
the branches of this tree, with a particular set of states of the 
nematic molecules. Each of the q branches of a Cayley tree 
with n layers of sites is built from q — 1 branches with n — 1 
layers. Due to the well-known thermodynamic pathologies 
associated with the surface of a Cayley tree 11211 . our aim is 
to study only those sites deep in the interior of the tree, which 
define the Bethe lattice. We thus focus on the properties of 
the central site (the root) of a tree with n surrounding layers as 
representative of the properties of the Bethe lattice obtained 
as n — > °°. The partition function of the MSZ model on such a 
tree is written as 



Z„ — Z n 



(4) 



where Z„ (p = x,y,z) is the partial partition function obtained 
when the molecule occupying the central site lies along the /j 
axis. 



The connections with the usual nematic structures of the 
liquid-crystalline systems are extracted from the tensor order 
parameter, 



l. 



J /JV, 



(5) 



where (■ • • ) indicates a thermal average, and the subscript 
refers to the central site. Terms of the form ((«o) 2 ) are related 
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Figure 1. Some layers of a Cayley tree with coordination q = 3. 
We also indicate some microscopic states of the sites of one of the 
branches of this tree. 



to the density of nematic molecules (in the Bethe lattice) with 
the symmetry axis along the /j direction (see Figure 1). We 
then adopt the correspondence 



<K) 2 > 

according to which we have 



lim 



Z>n 



(6) 



Un 2 + 2z„' 
satisfying the traceless property 

o^+Gf+e? =o. 

It is convenient to introduce the variables 
and 



(7) 



(8) 



(9) 



which are more adequate to distinguish between uniaxial and 
biaxial nematic structures. 

In order to calculate S„ and r|„, we notice that Z„ depends 
on the sum over states of the molecules in each of the q 
branches, 
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and 
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(10) 



(11) 



(12) 



where Q« is the partial partition function of a branch with n 
layers obtained when the molecule in the single site compris- 
ing the innermost layer lies along the fj axis. For these partial 
partition functions of branches, it is straightforward to write 
the recursion relations 



3 PA 



3 PA 



q-l 



(13) 



3 



Q„ 
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Defining the ratios 



p„ = and p£ = — , 



Eqs. O to ([151) lead to 



p5-i+pLi+ r 



9-1 



9-1 



(14) 



(15) 



= &(pJ!-i>p2-i) ( 16 ) 



with the derivatives calculated at the fixed point of interest. 

The trivial isotropic fixed point is stable for |Ai s0 | < 1, 
where 



Ai S o = (q-l) 



l-exp(-|PA) 



l+2exp(-|pA) 



(21) 



It is easy to find a numerical expression for the eigenvalue 
A un i, which is associated with the limit of linear stability of 
the uniaxial nematic fixed point, |A un ;| < 1. In Fig. [2] we 
draw graphs of A; so and A un ; as a function of temperature, for a 
typical value of the coordination q. The presence of a common 
range of stability requires a detailed thermodynamic analysis 
to choose the physical solution, corresponding to the smallest 
value of the free energy. The problem is further delicate, due 
to the need to avoid the pathologies produced by the surface 
sites of a Cayley tree. 



and 



Pn 



Pn-l+Pl-l + r 



9-1 



= «y(p5-i.tf-i). ( 17 > 



with r = exp(|PA). The connection between the variables 
(S„,r|„) and (p£,p«) is given by 



1 3 1 

5 = | 

2 2 [gM^n)}^+[gy{9 X n yn)}^+l 



11" = 



3 [gxfati)]^ ~[8y 

2 [&(pS»P£)]^+[ft>(pS,pi!)]^ 



, (18) 



(19) 
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The problem is now reduced to the analysis of the set of two 
nonlinear recursion relations ( fT6] l and ( fTTI i. Physical solutions 
on the Bethe lattice correspond to the (stable) fixed points of 
the mapping problem, which are solutions of 

p x = g x (9 X ,P y ) and p>'= ft .(p\p>). 

There is a trivial fixed point, p v = p v = 1 (S = T| = 0), cor- 
responding to a disordered (isotropic) phase, and an ordered 
fixed point, p x = p y ^ 1 (5^0 and r| = 0), corresponding to a 
uniaxial nematic structure. In this simple problem, it is easy to 
see that there is no possibility of a biaxial nematic fixed point, 
5^0 and r| ^ 0, which would require p* ^ p v . 

The occurrence of a first-order transition is associated with 
the existence of a range of temperatures in which both fixed 
points are dynamically stable (in other words, both fixed 
points can be reached from particular sets of initial condi- 
tions of the dynamic map IU5I1 ). This is indeed the case of 
the nematic-isotropic transition, as can be checked by a sim- 
ple linear-stability analysis, i.e. by finding the (degenerate) 
eigenvalues of the matrix 



/ 3 & 3 fo 
\ Tip 3p? 



(20) 




Figure 2. Temperature dependence of the eigenvalues associated 
with the analysis of linear stability of the fixed points of the simple 
MSZ model on a Cayley tree of coordination q = 3. The grey region 
indicates the joint range of stability (|A; so | < 1 and |A un ;| < 1). 



We now resort to a special technique ifTlllTrjll to find the free 
energy associated with the bulk attractors of the recursion re- 
lations. The idea of the method is to calculate that free energy 
by cleverly subtracting the contribution from the surface sites. 
According to an argument of Gujrati ifTTll . we first write the 
free energy of a Cayley tree with M layers, 



1 



■lnZ M 



P 

= i M) +qf[ M) + c l U l -\MV' r 
+q(q- 



f («) 



i) M -7j M) 



At-V 



(22) 



where Zm is the total partition function, defined in Eq. (O, 



AM) 



and /• is the free energy per site of molecules located at 
the y'th layer of the tree, while the coefficients of the /j M ' are 



4 



the number of sites in the j'th layer. We now rewrite this last 
equation as 



(M) 



M 



and notice that for a tree with M — 1 layers we can write 



M-l 



,(M-1) 



(«- i)f*_i = j +« e («- 

(The factor q—1 enters due to the fact that a M-layer tree has 
as many surfaces sites as q — 1 trees with M — 1 layers each.) 
As M — > oo, the free energies per site at the surface of both 
should become identical. In fact, we 



trees, f$*> and/^/' 
expect that f^-f^ 



for all values of j. Thus, 



F M -(q- 1)F M _, ^/,l M) +/( 



,(M-1) 



(23) 



Since /q M ' is the free energy associated with the central 
site of a M-layer tree, and the properties of this central site, 
for large M, are governed by the fixed points of Eqs. ( fT6l ) and 
( fTTI i. we should have 

jm/r^/r'w,, 

the free energy per site thus representing all sites in the 

Bethe lattice. Therefore, we conclude from Eqs. d22l and 
d23lthat 

fb = — tctt lim In — — r . 

2(3 (z m _^- 1 

We then use Eqs. (01 and dTOb-dTSb in the limit M 
the free energy 



. oo to write 



f b = f b (q,A/T;S,T]), 



(24) 



where S and r| are obtained from the fixed-point values of p r 
and p y via Eqs. cfT~8T > and (T% , and then determine which solu- 
tion corresponds to thermodynamic equilibrium. In Figure [3] 
we show a graph of the uniaxial order parameter, S, in terms 
of temperature. The dashed line corresponds to a coexistence 
of ordered (S ^ 0) and disordered (S — 0) solutions. 



III. LATTICE MODEL FOR A BINARY MIXTURE OF 
RODS AND DISKS 

Given a configuration {A,,} of rodlike (A,, = +1) and disk- 
like (A,, = — 1) molecular aggregates, the energy of the MSZ 
model for a binary mixture is written as 

#{3W} = -AE E hSfXjSf, (25) 

(i.j)^=x,y,z 

which leads to the canonical partition function 



z b = E' E ex P 

{h} {Iti} 



P^E E hhsfsf 

, (i,j)ti,V=x,y,z 



(26) 
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Figure 3. Graph of the order parameter 5 versus temperature (1/T = 
$Aq) for the MSZ model on a Bethe lattice of coordination q = 3. 
The dashed vertical line indicates a first-order transition. 



The sum over {A;} is restricted by the fixed concentrations of 
the molecular types, 



N r -N d = £A,,-, Nd =N — N r , 

i=\ 



(27) 



where N r (N c /) is the number of rodlike (disklike) molecules, 
and jV is the total number of molecules. It is now conve- 
nient to introduce a chemical potential and change to a grand- 
canonical ensemble, 

N 

E B = £ exp (pVV r ) Z b = £ £ ex P hPXff] , (28) 

N r =0 {X,} {%} 

with unrestricted sums over the sets of variables, and the ef- 
fective Hamiltonian 

X eS =-dt X -+ N )- A L E Wfsf. (29) 

1 V=i / (i,j)ftv=w 
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Figure 4. Two examples of configurations of rods and disks on the 
sites of a Cayley tree of coordination q = 3. 



Along the lines of the treatment of Section 2 for the simple 
MSZ model on a Cayley tree, we introduce some extra states 
to account for rod and disk variables (A, = ±1), and define a 



5 



set of six partial grand-canonical partition functions, Z„ x , E n " y , 
En' z , 2,i' x , En ,y , E«' z , associated with a tree of coordination 
q with n surrounding layers (see Figure 3J. The notation is 



such that E„' z represents the partial grand-canonical partition 
function of the tree with the central site in a state with A,, = 
+ 1 (rodlike) and n, along the z direction, n, ■ — ±(0,0, 1). In 
analogy to equations < TT~0b to dT2l . it is straightforward to write 



a r -3PA -3RA , 3BA , , 1 <7 

e-T + e^~ (q? + Q^) + Q** + (q** + q/' z ) , 



(30) 



3(3.4 
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3(5,1 



_t_ n d,x^ 



(31) 



(32) 
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(34) 
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+ e ^ (q/^ + Q^) + ^ c z + e*? (Q? + 
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(35) 



d v 

the quantities like Qj, ' now representing partial grand- 
canonical partition functions of a branch with n layers under 
the condition that the site in the innermost layer is occupied 
by a disk whose symmetry axis lies along the x axis. The full 
partition function of the tree is written as 



T?>Vt _i_ T7f,V _i_ 'h^Z _i_ ^d,x I ^d,y I ^d,z 



(36) 



We now use a similar parametrization as in Section 2. The 
tensor order parameter is given by 



3/1 
2\N 



N 



(37) 



where (■ ■ • } indicates a grand-canonical thermal average. The 
factor ((X, A,,) /AT) gives the difference between the concen- 
trations (number fractions) of rods and disks. It is then natural 
to introduce the correspondence 



— y ii \ -> iim — y 



(38) 



The ratio S« z /S n gives the concentration of rodlike particles 
with the symmetry axis along the z direction. Terms of the 
form ) 2 ) /N are recognized as the difference between 

the concentrations of rods and disks along the same fj direc- 



I 

tion, so that we have 
/ 1 



N 



P\2 



lim 



(39) 



As in Section 2, it is interesting to work with the tensor order 
parameter associated with the central site of a n-layer tree, 



O m = — y (Z ry - Z cLv ) + - 



2 ^ ' 



such that 



(40) 



(41) 



Again, it is convenient to introduce parameters S„ and r|„ to 
characterize the distinct phases. 

We can now write recursion relations for the branch partial 
partition functions Q^' v [formally obtained from Eqs. (f30b - 
© by taking E,?' v -> Q„ a - V , Q„ a ' v -»• Q^, and g -> g - 1], 
and rewrite them in terms of the set of five ratios 



d,x _ 0* 
Vn 



On' 



p d,y 
r n 



Qk 



</,v 



P " Q„ r ' : 



whose connection to the physical parameters S„ and T|„ is 
easily determined from Eqs. (f30b-(f35l). The resulting five- 
dimensional nonlinear mapping problem can be investigated 
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as in the previous section. The linear-stability analysis of the 
fixed points is then reduced to studying the eigenvalues of a 
5x5 matrix, analogous to that defined in Eq. (f2Qb . 

Depending on the values of temperature and chemical po- 
tential, there appears a biaxial nematic fixed point (5^0 and 
T| 7^ 0) besides two distinct uniaxial nematic fixed points (with 
5/0 and r| = 0); of course, there is also an isotropic fixed 
point (S — and r| = 0). In the fj x T plane (see Fig. |5), there 
is a large low-temperature region of stability of both uniax- 
ial fixed points. However, the fixed point associated with the 
biaxial nematic structure is dynamically unstable. 
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Figure 5. Stability lines in the /j x T plane for a mixture of rodlike 
and disklike molecules on a Cayley tree with coordination q = 3. 
Above the continuous line the isotropic phase is dynamically stable; 
below the dashed line at least one of the nematic uniaxial phases 
is dynamically stable; and between the dotted-dashed lines the two 
nematic uniaxial phases are dynamically stable. The biaxial nematic 
phase is dynamically unstable. 



The existence of common regions of stability of distinct 
fixed points requires the analysis of the associated free en- 
ergy in order to sort out the physically acceptable phases and 
to locate the coexistence (first-order) boundaries. Again, we 
use Gujrati's method to deal with the subtleties of a Cayley 
tree. As in Section 2, the free energy is written in terms of the 
function 

/r = fr(q,A/T^p r ^p^,p cL \p d ^p cU ). (42) 

The phase diagram for a tree of a typical coordination, q = 
3, is shown in Figure [6] At high temperatures, there is an 
isotropic phase (S = r| = 0). At low temperatures, there are 
two distinct uniaxial nematic phases, N r , with S > and r| = 0, 
with an excess of rods, and N c i, with S < and r| = 0, with an 
excess of disklike particles. In this thermalized formulation 
of the MSZ model, besides the biaxial nematic fixed point be 
dynamically unstable, it is associated with larger values of the 
free energy as compared to the uniaxial solutions, and cannot 
be thermodynamically acceptable for all coordinations q > 3. 
Of course, in the limit of infinite coordination (q — >• oo, A — >• 
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Figure 6. Phase diagram of the MSZ model for a mixture of rods 
and disks on the Bethe lattice (of coordination q = 3). We indicate 
two uniaxial nematic (N r and Nj) and an isotropic phase, which are 
bordered by first-order transition lines. 

0, with qA fixed), we recover all of the well-known results 
of the mean-field calculations. Finally, in Figure [7] we show 
the dynamic-stability lines of the isotropic and of the uniaxial 
nematic phases. As expected, the first-order transition line lies 
between the two stability lines. 
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Figure 7. Phase diagram presenting the nematic-isotropic first-order 
transition line and the two lines which indicate the dynamic-stability 
limit of these two phases. As expected the transition line lies in the 
region (shown in grey) where there is an overlap of dynamically sta- 
bility of the isotropic and the uniaxial nematic phases. 



From the experimental point of view, it is interesting to 
draw the corresponding temperature-concentration phase di- 
agram. We note that c = —dfi,/d/j, and perform a Legendre 
transformation to eliminate the chemical potential. In Figure 
|8] we draw the characteristic tie lines between coexisting uni- 
axial nematic phases. The regions of coexistence of isotropic 
and nematic phases are too narrow to be clearly represented in 
this phase diagram. 
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Figure 8. Temperature-concentration phase diagram of the MSZ 
model for a mixture of rods and disks on a Bethe lattice of coordina- 
tion q = 3. The tie lines indicate the coexistence of the uniaxial ne- 
matic phases, N r (rodlike molecules rich) and Nj (disklike molecules 
rich). There is no stable biaxial nematic phase. 



IV. CONCLUSIONS AND PERSPECTIVES 

We have formulated the problem of a Maier-Saupe- 
Zwanzig model on a Cayley tree as a set of nonlinear discrete 
recursion relations, whose attractors correspond to physically 
acceptable solutions on the Bethe lattice (deep in the interior 
of the tree). Due to the presence of first-order transitions, we 
have large regions of overlap of stability of distinct attractors. 
We then resort to an ingenious scheme to obtain the Bethe- 
lattice free energy associated with the coexisting attractors and 
to determine the thermodynamically acceptable solutions. 

We first considered the simple MSZ model on a Cayley tree. 
In this problem, there are just two fixed points, corresponding 
to disordered and uniaxial ordered structures, which are both 
(dynamically) stable in a certain intermediate range of tem- 
peratures, and which is an indication of occurrence of a first- 
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order transition. We then use Gujrati's method to locate the 
first-order boundary between the high temperature isotropic 
and the low temperature uniaxial nematic phases. The MSZ 
model for a binary mixture comes from the introduction of a 
new set of shape variables and the definition of an effective 
Hamiltonian in the grand-canonical ensemble. The problem 
is more involved, but we can perform a relatively simple sta- 
bility analysis of the fixed points, and draw a phase diagram in 
terms of temperature and chemical potential. At low temper- 
atures, we find large regions of overlap of stability of fixed 
points associated with uniaxial phases. However the biax- 
ial structure is unstable both from a dynamic-stability and a 
thermodynamic standpoint. The main qualitative features of 
the phase diagrams on the Bethe lattice agree with previous 
mean-field predictions ||8|] . 

The framework employed in the study of the binary mixture 
can be extended in straightforward ways to deal with more 
complicated interaction potentials, although of course requir- 
ing greater analytical effort. An obvious extension is to con- 
sider a pair interaction described by 

£[U; + &(^+^) + c]Sfsf, 

which breaks the symmetry between rod-rod and disk-disk 
interactions. It is also possible to introduce dilution and 
isotropic repulsive forces to mimic the interaction potentials 
leading to stable biaxial phases in recent computer simulations 
|0]. We hope to report on such extensions in future publica- 
tions. Other extensions possibly leading to stabilization of a 
biaxial pha se in mixtures, such as the introduction of polidis- 
persity il7lll8ll . can also be considered, at least in principle, 
although the analytical work becomes quite involved. 
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